Solidification in a channel 
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We simulate solidification in a narrow channel through the use of a phase-field model with an 
adaptive grid. In different regimes, we find that the solid can grow in fingerlike steady-state shapes, 
or become unstable, exhibiting unsteady growth. At low melt undercoolings, we find good agreement 
between our results, theoretical predictions, and experiment. For high undercoolings, we report 
evidence for a new stable steady-state finger shape which exists in experimentally accessible ranges 
for typical materials. 
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The classic example of pattern selection in driven sys- 
tems is the shape of a dendritic crystal's tip, which is 
growing into a supercooled melt, mediated by the diffu- 
sion of the latent heat of solidification p]J^| . Simple theo- 
ries, which neglect the effect of surface tension k, predict 
that only the product of the velocity V and dendrite tip 
radius p is determined for a given undercooling. How- 
ever, the modern theory of microscopic solvability, which 
incorporates the surface tension, predicts the unique op- 
erating condition for a steady-state tip, provided that 
k is anisotropic 0. This work generalized earlier work 
on local models of solidification ||^|] and, in particular, 
work on the closely related problem of viscous fingering 
in Hele-Shaw cells M . The latter phenomena corresponds 
to the fingering instability which occurs when a relatively 
inviscid fluid, such as water, is forced by a pressure gra- 
dient to invade a more viscous fluid, such as oil. It can be 
shown that, if this process takes place between two glass 
plates which form a long channel of width L — called 
a Hele-Shaw cell — a long steady-state finger forms and 
grows §. As in the case for dendritic tips, in the absence 
of surface tension, theory predicts fingers of any width. 
The theory of microscopic solvability predicts a unique 
finger width, e.g., in units of the channel width, the fin- 
ger width approaches A = 1/2 in the limit of k — > 0. 
Mathematically, both the dendrite tip and the viscous 
finger are described by similar equations, namely, diffu- 
sion or the Laplace problem in the bulk phases, subject 
to the boundary conditions of the moving interface, and 
the driving force at infinity far from the interface. In the 
limit in which the diffusion length I u becomes very large 
compared to the channel width, where Id — D/V and D 
is the diffusion coefficient (that is, if the Peclet number 
P = L/ln is small), the two problems become almost 
identical. This close relationship has led to many beau- 
tiful insights into both phenomena For example, 
the "free" viscous finger generically shows unsteady tip 
splitting behavior, which can be attributed to the fluid's 



isotropic surface tension. However, creating an effective 
anisotropic surface tension, by scratching one of the glass 
plates, stabilizes the tip, giving rise to dendritic patterns 

Here we report the results of a theoretical investigation 
into the alternate case: solidification constrained to grow 
in two dimensions within a long channel. For small P, 
this is closely related to viscous fingering with anisotropic 
surface tension. We compare our numerical results for a 
phase-field model of solidification, with no adjustable pa- 
rameters, to theory Q and experiment finding good 
agreement, thereby providing the first quantitative test 
of solvability theory in this context. For larger P, the 
conservation law for the latent heat plays an important 
role. Indeed, for large undercoolings, we find evidence for 
a new steady-state finger, which we relate to local models 
HQ of dendritic growth. 

Our approach takes advantage of recent advances in 
computer power, as well as algorithmic development. We 
use a two-dimensional phase-field model of solidification, 
|lo| fl8[ where a scalar field <j)(x, t) determines whether 
the phase at field point x and time t is solid {4> = +1) or 
liquid (cf> = — 1). A useful approach introduced in Ref. 
|]l5f allows for the elimination of kinetic undercooling. 
As well, we use finite-element adaptive-grid refinement 
techniques, introduced by one of us and collaborators 
|]l6f , wherein a fine mesh is used near the rapidly varying 
solid-melt interface, and a coarse mesh is used in the bulk 
phases where quantities vary slowly with space and time. 
Let the dimensionless temperature be = c p (T — Tm) / 
where c p is the specific heat at constant pressure, Tm is 
the melting temperature, and £ is the latent heat of the 
fusion. Then the equations of motion are |l6| , |ri| 

and, 
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where D = Vtq/Wq, T> is the thermal diffusivity. All 
quantities were measured in units of the microscopic 
time and length scales, tq and Wq. The parameter 
7 is chosen to be 1.5957D to eliminate kinetic under- 
cooling [fL5|| , and thereby relate To and Wo to experi- 
mental values. Anisotropy enters through A(n), where 
n = V</>/|V0|, is the unit vector normal to the con- 
tours of 0: A(n) = [1 - 3e][l + ^((n,) 4 + Kp], 
where e is a constant [fuj[ . We used e = 0.05 |p|,|l9|. 
The width of the interface is W(n) — WoA(n) and the 
characteristic time r(n) = ToA 2 (n) [|l5). Using d a = 
0.5539/D, the boundary condition at the interface fl§| l 
is ®int = —d(n)K, where K is the interface curvature, 
and d(n) — do[A(n) + d 2 A/d(cos~ 1 n x ) 2 ] is the capillary 
length. The simulations were performed using Neumann 
boundary conditions: n ■ V<^> = and n ■ V0 = 0, so that 
heat cannot leak out of the channel. Solidification is ini- 
tiated by a small half disk of radius Ro centered at the 
bottom of the channel. The interface is initially set to its 
equilibrium value <p(x) = — tanh(x — Ro)/\/2. The initial 
temperature is = in the solid phase, which decays 
exponentially to its asymptotic value = A at x — > oo 
in the far field. The time step was chosen to be 0.016, 
and the minimum grid size for the adaptive grid was 0.78 
and 0.39, for the low and high P regimes, respectively. 

With this approach, we can make a quantitative com- 
parison to experiment as well as theory, since as P — > 
these equations reduce to those of the anisotropic vis- 
cous finger. Motivated by the experiment of Molho, Si- 
mon, and Libchaber p], the first set of simulations were 
done in a wide channel L/do — 121, 000 with rectangular 
cross section of 1600W-0 by 102400Wo (319.3 /imx2 cm), 
similar to the experimental capillary cells || (200/im x 6 
cm). To examine the low Peclet regime, the simulations 
were done at low undercoolings A = 0.01, 0.05, 0.1, 0.125, 
0.15, with D = 40 and Rq = 30. After a short period 
of acceleration, tip velocity decreases to an asymptotic 
value, following approximately an inverse square root of 
time. To estimate this asymptotic tip velocity and width 
of the fingers, we extrapolated our transient results us- 
ing a scaling ansatz, as in Ref. and by fitting to a 
high-order polynomial. The scaling ansatz is shown as 
an inset in Fig. 1; we scaled the width of the finger with 
t a , and the length of the finger with t 13 , obtaining very 
good data collapse for the entire shape of the finger for 
different undercoolings p(| . 

This method and the polynomial fits gave consistent 
results for A, which are shown in Fig. 1, where V appears 
through the dimensionless surface tension n — gswQ^p . 
We found that for these channel dimensions the dimen- 



sionless finger width, A, is smaller than 0.5, and decreases 
with increasing velocity, as predicted in the solvability 
theory of Kessler, Koplik, and Levine Q]. Also shown 
are the results of the experiment, and the solvability solu- 
tion for the anisotropic viscous finger Qj2l|,|2^] . With no 
adjustable parameters, our phase-field model results are 
consistent with the experimental observations ||pl|and 
compare very well to the theoretical prediction ]7|, |2l| , p2[ . 
Our results provide the first quantitative test of that pre- 
diction of solvability theory for anisotropic viscous fin- 
gers, confirming those predictions. 

This correspondence between solidification in a chan- 
nel and the anisotropic viscous finger, which we have ex- 
ploited above, is restricted to small Peclet numbers, and 
to a transient time regime which diverges as P — » 0. 
For later times, and larger P, the situation is different: 
from thermodynamics, the only asymptotic solution for 
an arbitrarily long solid finger within a liquid channel 
is A = A, since A is the fractional amount of solid in 
equilibrium @|||. For long times, far beyond those in- 
vestigated above, it is expected that the finger shape ob- 
tained above is unstable due to latent heat piling up in 
front of the finger. For example, Pelce [^3) has pointed 
out that the finger is unstable to a tip widening instabil- 
ity unless dX/dV > 0: If A is a decreasing function of V, 
then a perturbation giving rise to a wider finger causes 
the finger to slow down, leaving more heat in the vicinity, 
causing it to slow down more, and flatten out until the 
finger reaches the size of the channel. 

To investigate this regime, we simulated larger Peclet 
numbers in a narrow channel L/do = 355, of dimension 
100Wo by 1600Wo- The simulations were performed at 
A = 0.25 - 0.95, with D = 2 and R = 4. Fig. 2 shows 
the dependence of V on time. There is a marked quali- 
tative change at A w 0.5. For A < 0.5, a tip widening 
instability takes place. If we estimate an effective V and 
A from a transient regime before this instability takes 
place, we find that A is a decreasing function of V, as 
would be expected: the dependence of A on undercooling 
for these fingers is shown by the open circles in Fig. 3; 
inset (a) of Fig. 3 shows one of these fingers while the 
instability is taking place. 

In contrast, for A > 0.5, in the presence of anisotropic 
surface tension, the fingers rapidly achieve a well-defined 
constant velocity (closed circles in Fig. 2) and shape. Fits 
to V give P = 5.4 — 49. In Fig. 3, we show that their 
finger widths obey A = A; inset (b) shows one of these 
fingers. Hence these fingers satisfy thermodynamic and 
stability requirements, since A cx V. This new class of 
stable fingers fingers becomes unstable to tip splitting 
if the surface tension is isotropic. We do not have an 
explanation for the apparent stability criterion A > 0.5, 
although it is interesting to note the reappearance of this 
bound for isotropic viscous fingers in the new context of 
rapidly growing anisotropic solid fingers. Nevertheless, 
for high undercoolings close to unity, we can suggest an 
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origin for these new fingers. 

For high undercoolings, the thermal boundary layer 
around the solidifying finger becomes small. As shown 
by Ben- Jacob, Goldenfeld, Langer, and Schon in 
the limit A — > 1 the interface velocity and curvature 
are determined locally, as in, e.g., the geometric || and 
boundary-layer 0] models of solidification. In such mod- 
els, the steady-state shape is determined by the curva- 
ture K as a function of only the normal velocity in the 
limit of vanishing surface tension. That is, for small K 
and n x , K — f(n x ) — an x + bn^, + • ••. For the free 
boundary-layer model |3j , the coefficient a is exactly zero, 
and b oc (1 — A), giving rise to a parabolic shape. How- 
ever, it is easy to show that a nonzero a is required for 
a finger of finite width; here the value of a is determined 
by the thermodynamic requirement A = A. This gives 
rise to the new steady-state solution: 

x = -B(A/n) lncos(7n//A) 

where 1/B = y/l + (6A/2tt) 2 - (bA/2%). This form gives 
good one-parameter fits to the fingers as shown in Fig. 4 
(and fits which are significantly better than those from 
the Saffman- Taylor || shape 6 = 0). Since the form is 
obtained from an expansion in n x , it is a better fit in the 
tails that at the tips. Nevertheless, the fitted dependence 
of b on A is reasonably well described by b oc (1 — A), for 
A — > 1 as expected. 

To conclude, we have investigated solidification in a 
channel in low and high undercooling regimes. Our re- 
sults for low undercooling are consistent with earlier ex- 
perimental work, and compare well to theory, provid- 
ing the first quantitative test of solvability theory for 
anisotropic viscous fingers. For large undercoolings, we 
find new phenomena, a solid finger which satisfies stabil- 
ity and thermodynamic criterion. We further argue for 
an analytic form of the shape, based on local models of 
solidification, which fits our results from numerical sim- 
ulation. We note that this new "finger" can be studied 
further experimentally and theoretically. 
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FIG. 1. Finger width A vs. velocity V for small Peclet 
numbers. Solvability theory [7,21,22] (solid line, e = 0.05), 
experimental results [8] (empty circles), and present simula- 
tions (filled circles, e = 0.05). Inset: scaling collapse used 
to estimate V for A = 0.05 (thin line), and A = 0.15 (thick 
lines). X s is scaled length of finger, Y s is scaled width. For 
A = 0.05 all profiles collapse. For A = 0.15, seven profiles 
are shown, which collapse quickly for late times. 

FIG. 2. Velocity V vs. time for various undercoolings 
A, for larger Peclet numbers. Fingers with tip- widening in- 
stability A < 0.5 (open circles). Note qualitative change in 
velocities for stable fingers A > 0.5 (filled circles). 

FIG. 3. Finger width A vs. undercooling A for larger 
Peclet numbers. Fingers with tip-widening instability (open 
circles, inset (a)). Stable fingers (filled circles, inset (b)). 

FIG. 4. Numerical (open circles), and analytic (lines) re- 
sults for stable finger tip profiles, at undercoolings A = 0.55 - 
0.95. Analytic fits, from an expansion valid in the tail, involve 
one parameter which should obey b oc (1 — A) for A — > 1. b 
vs. A (bottom inset) shows reasonable agreement. 
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